MIT-CTP 4108 



Quantizing Majorana Fermions in a Superconductor 

C. Chamon 1 , R. Jackiw 2 , Y. Nishida 2 , S.-Y. Pi 1 , L. Santos 3 

1 Physics Department, Boston University, Boston MA 02215 
2 Physics Department, MIT, Cambridge MA 02139 
3 Physics Department, Harvard University, Cambridge MA 02138 



Abstract 

A Dirac-type matrix equation governs surface excitations in a topological insulator 
in contact with an s-wave superconductor. The order parameter can be homogenous or 
vortex valued. In the homogenous case a winding number can be defined whose non- 
vanishing value signals topological effects. A vortex leads to a static, isolated, zero 
energy solution. Its mode function is real, and has been called "Majorana." Here we 
demonstrate that the reality/Majorana feature is not confined to the zero energy mode, 
but characterizes the full quantum field. In a four-component description a change of 
basis for the relevant matrices renders the Hamiltonian imaginary and the full, space- 
time dependent field is real, as is the case for the relativistic Majorana equation in the 
Majorana matrix representation. More broadly, we show that the Majorana quantiza- 
tion procedure is generic to superconductors, with or without the Dirac structure, and 
follows from the constraints of fermionic statistics on the symmetries of Bogoliubov-de 
Gcnncs Hamiltonians. The Hamiltonian can always be brought to an imaginary form, 
leading to equations of motion that are real with quantized real field solutions. Also 
we examine the Fock space realization of the zero mode algebra for the Dirac-type 
systems. We show that a two-dimensional representation is natural, in which fermion 
parity is preserved. 

Introduction 

Majorana bound states arise as zero energy states in two-dimensional systems involving 
superconductors in the presence of vortices. [1, 2, 3, 4] These zero modes have attracted 
much attention recently, in part because of the possibility that they can realize "half- 
qubits" within topological quantum computing schemes [3]. The basic idea is that two far 
away Majorana bound states, real fermions, can be put together into a complex fermion 
acting on a two-dimensional Hilbert space spanned by the states |0) and |1). Hence, two 
Majorana fermions comprise one qubit, which is protected against the environment if the 
vortices binding the Majorana fermions are kept far away from each other. 
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The first example of a zero mode in a two-dimensional superconductor was presented in 
Ref. [1]. More recently it has been stated that the proximity effect at the interface between 
an s-wave superconductor and the surface of a topological insulator can be described by a 
planar Dirac equation [4] , providing a physical realization of the mathematical structure of 
Ref. [1] Other examples of Majorana bound states arise in systems with a non-relativistic 
kinetic term and a p± = p x ± ip y interaction with a vortex order parameter (i.e., p-wave 
superconductors). [2, 3] These types of bound states have been the subject of much recent 
interest [5, 6, 7, 8]. The focus of the discussions of Majorana fermions in superconductors 
have focused thus far on the zero modes. 

However, Majorana's original work [9] was actually quite more general, and did not 
address a single mode but instead a whole field. What he showed was that it was possi- 
ble to construct a representation of the Dirac equation that admits purely real solutions. 
The particles that follow from his construction are their own anti-particles, and thus nec- 
essarily neutral. What was striking about Majorana's proposal was that these particles 
were fermions - bosonic neutral particles represented by real fields are common, pions and 
vector bosons, such as photons, being simple examples (see Ref. [10] for a perspective on 
Majorana fermions). 

In this paper we look at three issues regarding the quantization of Majorana fermions, 
beyond simply the zero modes, in superconductors. First, we look specifically at the 
case of Dirac-type systems describing s-wave induced superconductivity on the surface of 
topological insulators. There, we find that the entire ip field of the superconductor model 
(and not merely particular modes) obeys equations that are analogous to the Majorana 
equations of particle physics. The equations of motion for the fields can be brought to a real 
form, and the fermionic solutions are real and therefore their own anti-particles. Indeed, 
other than the fact that surface states are 2D, the topological insulator-superconductor 
system can be brought to the exactly same form that was discussed in Majorana's original 
formulation of real relativistic fermions. 

Second, we note various topological features of the Dirac-type model. We compute the 
Pontryagin index associated with the k-space dispersion, and find it to be ±1/2, which is 
an indication of the existance of zero modes in the presence of vortices. We then present the 
Fock space level structures that accommodate an isolated, zero energy state, which arises 
in the presence of a vortex. In particular we show that fermion parity can be preserved, 
even with a single zero energy state. As we discussed above the Majorana zero modes 
are usually thought of "half" qubits, as two of them make up a complex fermion with a 
two-dimensional Hilbert space. Here we ruffle this simple view by quantizing the theory 
in the infinite plane in the presence of a single vortex. A sole Majorana zero mode exists, 
but a two-dimensional Hilbert space remains. In a finite system, another zero mode would 
appear at the edge, which is however absent in the infinite plane. 

Third, we show that the Majorana quantization procedure that we discuss for the Dirac- 
type equations describing s-wave induced superconductivity on the surface of topological 
insulators does extend, more broadly, to any superconductor. A description of Bogoliubov- 
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de Gennes Hamiltonians using Majorana modes has been noted by Senthil and Fisher [11] 
for systems where spin rotational symmetry is broken (classes D and Dili of Ref. [12]). 
Here we show rather generically that lack of spin rotation symmetry is not a necessity, and 
thus classes C and CI of Ref. [12]) also realize Majorana fermions. All one actually needs 
is to have fermions, and hence these results hold for any superconducting system made 
of half-integer spin particles, regardless of the size of the spin. What we show is that the 
constraints imposed by fermionic statistics on the symmetries of Bogoliubov-de Gennes 
Hamiltonians always allow one to bring the Hamiltonian in the Nambu representation to 
an imaginary form. In turn, Schrbdinger's equation with this imaginary Hamiltonian leads 
to a real equation of motion for the fields, as in Majorana's construction. The real field 
solutions in the constrained doubled Nambu space can then be quantized as Majorana 
fields. 
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Quantum Structure of the Superconducting Model 



Let us start by analyzing the planar Dirac-type systems realized on the surface of a topolog- 
ical insulator, placed in proximity to an s-wave superconductor. The Hamiltonian density 
for the model under discussion acts on two spatial dimensions. [4, 1] 



H = ij)*p- fa + fa p + fa - p(ip* fa + fa fa) 

+ a fa fa + A*^ t (i) 

Here fa tl are electron field amplitudes, p± = —id x ± d y , p, is the chemical potential (which 
was omitted in the ref. [1]) and A(r) is the order parameter that is constant in the homoge- 
nous case or takes a vortex profile in the topologically interesting case: A(r) = v(r) e %9 , in 
circular coordinates. Equivalently, in a two-component notation 

il = </>* (<r • p - fa + - A fa io% ifi* - - A* fa ia% fa . (2) 

Now if) = f ^ ^ , and <x comprises the two Pauli matrices (a 1 , a 2 ). The (2+l)-dimensional 
equations of motion for (1), (2) 

id t fa = p- fa - fi fa + A fa ^ 
id t t/) i = p+ V't - M ipi ~ A 

can be presented in two-component matrix notation. 

i d t ip = (<r • p - /x) tp + A icr 2 ?/>* (4) 

When the chemical potential is absent, and A is constant, the above system is a (2+1)- 
dimensional version of the (3+l)-dimensional, two component Majorana equation, which 
in (3+l)-dimensional space-time describes chargeless spin 1/2 fermions with "Majorana 
mass" | A |. [13] 

A static solution to (3), equivalently (4), with a vortex profile for A, can be readily 
found. It corresponds to a zero energy mode. With / and g real in the Ansatz 

■0t = /(r) exp {— i 7r/4 — y(r)} 

fa = g(r)exp{i(6 + ir/4:)-V(r)} (5) 
V'(r) = v(r) 



(3) reduces to 



(r 9 )' = »rf 
f = -y-g ■ 
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(Dash signifies r - differentiation). Regular solutions are Bessel functions 



f(r) 

a fa) 



NJo far) 
N,h far) 



(7) 



with ./V a real normalization constant. [14] 

While the static, zero energy mode is readily obtained from eq. (3), for the finite energy 
modes, we must take account of the fact that ip^^ip^ mix with their complex conjugates. 
Therefore, one cannot separate the time dependence with an energy phase. Correspond- 
ingly one cannot construct a Hamiltonian energy eigenvalue problem, which is the usual 
first step in the quantization procedure. 

Progress is achieved by doublings the system with a four-component spinor. 



V -V> t * J 



ia 2 V>* 



(8) 



An extended Hamiltonian density % leads to equations for which are just two copies of 
(3) or (4). 



H 





1 



0) 



Here T denotes transposition. Because the last two components of ^ are constrained by 
their relation to the first two, ^ satisfies the (pseudo) reality constraint 



(10) 



with C = C~ 



C* = C 1 



i a 



-i a 




To proceed, one ignores the constraint (10) on ^, and works with an unconstrained 



four-spinor $ 



. Time can now be separated with the usual phase Ansatz, and the 



energy eigenvalue spectrum can be found 

h® E 



id t <5> , $ = e- iEt $ E 



E<S> E 



(11) 



These are the Bogoliubov-de Gennes equations for the superconductor problem. In the 
particle physics application, the unconstrained four-component equation is just the Dirac 
equation describing charged spin 1/2 fermions. When the (pseudo) reality constraint is 
imposed, one is dealing with the four-component version of the Majorana equation. [13] 
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Observe that h in (9) possesses the conjugation symmetry 

C~ l hC = -h*, (12) 

which has the consequence that to each positive energy eigen mode there corresponds a 
negative energy mode. 

C$\e = ®-e (13) 

A quantum field may now be constructed by superposing the energy eigen modes §e 
with appropriate creation and annihilation quantum operators. It is here that we again 
encounter the Majorana construction: the unconstrained fermion four-spinor $> is like a 
"Dirac" fermion spinor, governed by a Hamiltonian, which satisfies a conjugation symmetry 
(12) that leads to (13). Then the spinor which satisfies the (pseudo) reality constraint 
(10), is like a "Majorana" spinor, viz. a "Dirac" spinor obeying a (pseudo) reality condition. 

With the eigen modes one can construct a quantum field <I>. It can be an unconstrained 
"Dirac" field operator. 

-iEt (Tl „ , V 1 l,t „-iEt , 



E<0 

~ iEt $e + J2 " iEt 
E>0 E>0 



E>0 E<0 

E a E e- iEt <S> E + E b E e iEt C$ E 1 ' 



Here the a E operator annihilates positive energy excitations (conduction band) and the 
b E operator creates negative energy excitations (valence band). Since 6 is unconstrained, 
a and b are independent operators. Their conventional anti-commutators ensure that the 
unconstrained fields satisfy Dirac anti-commutation relations. 

{^(r),^!-')} = (15a) 
$t(r')} = <%£(r-r') (15b) 

For the superconductor /topological insulator system under consideration $ — > fy, and 
the quantum field ^ satisfies the constraint 

Cij *t = %. (16) 
This is achieved by setting b = a in (14). 

(a E e- iEt <S> E + a E e iEt C<S> E ^J (17) 

Owing to the constraint (16) the anti-commutators take a "Majorana" form. 

{^(r),^(r')} = a^(r-rO (18a) 
{^(r),^(r')} = M(r-r') (18b) 
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These also follow from (17), with a^,a^ obeying conventional anti-commutators. We have 
ignored possible zero-energy states; they will be discussed at length below. 

In the final result (17), i& retains the Majorana feature of describing excitations that 
carry no charge. This is true for the entire quantum field not only for its zero energy 
modes (if any), which are emphasized in the condensed matter literature. Explicitly we 
see this by examining the conserved current that is constructed with the unconstrained 
"Dirac" field $. 





cr 





4>; 



(19) 



When the above is evaluated on the constrained field all terms vanish. This is to be 
expected for a Majorana field which carries no charge. 

One may also consider a chiral current constructed with the "Dirac" field 



(p5, Je 



But with non-vanishing A this is not conserved. 







' I 




a " 




V 


-I 




cr 








ij 




V / 



(20) 



Of 



p 5 + V • J 5 



-2i 



A 
-A* 



These results persist when the constraint (16) is imposed on — > 



0_ 

Ot 



(21) 

(22) 
(23) 



Thus no conserved current is present in the superconductor model (1). 

The Major ana/reality properties are obscured by the representation of the Dirac ma- 
trices employed in presenting the 4x4 Hamiltonian h (9). As written, the matrices in h 
are given in the Weyl representation. 



a 



cr 
-cr 



P 



/ 
/ 



75 



/ 
-/ 



h = ol ■ p - n 75 + j3 A R - if} 75 A/ 



(24) 
(25) 



A^/ are the real and imaginary parts of the order parameter. One may pass to the 
Majorana representation by conjugating with the unitary matrix 



V 



Q- Q+ 
Q+ -Q- 



Q ± = \ {l±a 2 ) 



(26) 
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Then h becomes 

v -i hV _f-Py Ps^ + iAjN //id 2 -A R a 2 

yPz^-iAj p y J \ -A R a 2 -ilo 1 

This is manifestly imaginary, and the conjugation matrix C in (10) becomes the identity, 
so that the (pseudo) reality constraint on \I/ becomes a reality condition. [17] 
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Homogenous Order Parameter 



For constant A = me tuJ , we pass to momentum space with an e* kr Ansatz in (11). The 
energy eigenvalue is 

E = ±^J{k ± fi) 2 + m 2 (28) 

with no correlation among the signs. For fixed k there are two (±/i) positive energy 
solutions and two negative energy solutions. They become doubly degenerate at /i = 0. 

e iw a 3 

The degeneracy occurs because at /i = 0, h commutes with S 



e~ iuJ a 3 

when the phase a; of A is constant. The energy (28) is non- vanishing for all values of the 
parameters; there is no zero-energy state. 

The operators a^,a^ and the eigen modes 3>e, which are explicitly presented in Ap- 
pendix A, are labeled by the momentum k and a further (+, — ) variety describing the two- 
fold dependence on fj, of E± = \/(k^f fi) 2 + m 2 ^ \m\. The quantum operator \f is con- 
structed as in (17), which with notational changes [cle — > a n (k); $b — > <5 n (k); n = (+, — )] 
reads explicitly 

*(*.*) = E/(0|«n(k)e-^ i - k - r )$ n (k) 

" L (29) 

with positive energy eigenfunction $ n (k) carrying energy eigenvalue E n . The conjugation 
condition (13) now states that C<fr*(k) is a negative energy solution at (— k). 

Actually we can suppress the lower two components of $ n (k) in (29), because they 
repeat the information contained in the upper two components, owing to the subsidiary 
condition (16). In this way from the four-component spinors recored in Appendix A, we 
arrive at a mode expansion for the electron field operators ■ 



d 2 k ( 1 / ,. , v.l A 



a+(k) e -W+t--kT) 



(2vr)2 \^ET\ y/E + -k + f , ' 

- a f + (k) e'^+^-k-r) -k + n e-*\-) J 

+ ^L= fa_(k) e -^ E -*- k ' r ) 



V y/El + k + (1 

+ al (k) e ''^-*- k ' r) V^- + fc + M |+) ) } 



(30) 
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Here |±) are the 2-component eigenvectors of cr ■ k : |+) = 




where tp arises as k = (cos ip , sin (p) . One verifies that (30) satisfies (3). 

The "Majorana" character of this expression manifests itself in that the particle anni- 
hilation operators a±, associated with the positive energy eigenvalues E±, are partnered 
with their Hermitian adjoint creation operators a^_, which are associated with the negative 
energy — E± modes. By contrast, for a "Dirac" field the negative energy modes are asso- 
ciated with the anti-particle creation operators which anti-commute with a±,a±. In 
other words, in the Majorana field operator (30) the anti-particle (hole) states are identified 
with the particle states. 
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Topological numbers 



When n is absent and S commutes with h, we may equivalently work with b! = S h, which 
possesses the same eigenvectors as h, common with the eigenvectors of S. However, h' has 
the appealing form 

ti = T, a n a (a = 1,2,3). (31) 

Here n l = k l (i = 1, 2) and n 3 is A e~ %w = m, i.e the constant phase of A is removed, so m 
is a real constant, but of indefinite sign. The matrices E a 

( -ie*"e«V \ _ / a 3 0\ 

satisfy the SU(2) algebra, as is explicitly recognized after a further unitary transformation. 

tr's.*-^") (33) 



^[ e -4 p + " J ' p± ^ 2 (/±(J3) (34) 

When the Hamiltonian is of the form (31), we can consider the topological current in 
momentum space. [20] 

Kil = ^ e abc h a d a h b dph c (n = n/|n|), (35) 
and evaluate the topological number by integrating over 2-dimensional k-space. 

M = [ d 2 kK°(k) = — [ d 2 k - - ^— — = (36) 
J ^ ' 87T J (fc2 +m 2)3/2 2 |m| v ' 

The non- vanishing answer ±1/2, depending only on the sign of m, is evidence that the 
model belongs to a topologically non-trivial class. It is also a hint that topologically 
protected zero modes exist in the presence of a vortex. [Although vortex based zero modes 
are also present for fi ^ 0, we do not know how to define a winding number in that case.] 
We can understand the fractional value for N . The unit vector 



h a = (k cos (p, k sin ip, m) /y k 2 + m 2 (37) 

maps gins at k = 0, n a is at the north or south pole 

(0,0, ±1). As k ranges to oo, h a covers a hemisphere (upper or lower) and ends at the 
equator of S^ 2 \ Thus only one half of is covered. 
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Single Vortex Order Parameter 



When A takes the vortex form, A(r) = v{r)e td , eq. (11) possesses an isolated zero energy 
mode 

ia 2 V>o 

with t/?Q determined by (5) and (7). Note that 



*o=( . I" (38) 



C%* = % (39) 

There are also continuum modes. 

The operator field \l/ is now given by an expansion like (17), except there is an additional 
contribution due to the zero mode controlled by the operator A. 

* = E {a E e- iEt <S>E + a ] E e iEt C<S>* E ) + Ay/2% (40) 

E>0 

(The y/2 factor will be explained later.) Due to (16) and (39), A is Hermitian A = A\ 
anti-commutes with (aE,a E ) an d obeys 

{A,A} = 2A 2 = 1. (41) 

The question arises: How is A realized on states? Two possibilities present themselves: 
two disconnected one-dimensional representations, or one two-dimensional representation. 

In the first instance, we take the ground state to be an eigenstate of A. The possible 
eigenvalues are ± so there are two ground states, |0+) with eigenvalue +^75, and |0— ) 

with eigenvalue — 4j. No local operator connects the two, and the two towers of states 
built upon them 

4 4/ a L---l 0± ) 

define two disconnected spaces of states. Moreover, one observes that A has a non- vanishing 
expectation value (0±| ^4 |0±) = Since A is a fermionic operator, fermion parity is 

lost. [21] 

In the second possibility, with a two-dimensional realization, we suppose that the vac- 
uum is doubly degenerate: call one "bosonic" the other "fermionic" |/), and A connects 
the two 

A\f) = ± \b) (42) 

A\b) = ± I/) . (43) 
(Phase choice does not loose generality.) 
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Again there are two towers of states 



E a Ei a En ■ ■ ■ I/) ' a E a El a Ell 



...\b) 



but now A connects them. With this realization, fermion parity is preserved when \b) 
and |/) are taken with opposite fermion parity. Of course since A is Hermitian, it can be 
diagonalized by the eigenstates. 



This regains the two states of the two one-dimensional realizations. But the combination 
\b) ± |/) violates fermion parity as it superposes states with opposite fermion parity. 

There does not seem to be a mathematical way to choose between the two possibilities. 
But physical arguments favor the fermion parity preserving realization. First of all, there 
is no reason to abandon fermion parity; if possible it should be preserved since it is a 
feature of the action. Also arguments against combining states of opposite fermion parity 
may be given: Since bosons and fermions transform differently under 2ir- spatial rotations; 
the |0±) states in (44) are not rotationally covariant, but transform into each other. [This 
argument is completely convincing in a (3+l)-dimensional theory. In (2+1) dimensions 
the anyon possibility clouds the picture, and in (1+1) dimensions the argument cannot be 
made, because spatial rotations do not occur.] Furthermore, time inversion transformations 
work differently on bosons and fermions: T 2 is / for bosons and —I for spinning fermions. 
The superposed states (44) are not invariant under T 2 , rather they transform into each 
other. [This argument can be made for (2+1) dimensional models, but in (1+1) dimensions 
spin is absent so the fermion parity violating option cannot be ruled out. Furthermore, 
the fermion-boson equivalence of (l+l)-dimensional models obscures the status of fermion- 
boson mixing. Indeed it is argued within super-symmetry that fermion parity is lost in the 
presence of solitons in (1+1) dimensions "due to boundary effects." [22]] 

[Any argument based on time inversion transformations requires viewing the complex 
valued vortex configuration as arising from the degrees of freedom of an enlarged model, in 
which the vortex emerges from the dynamics of the extended model (Abrikosov, Ginzburg, 
Landau). Otherwise, a vortex background is not T-invariant.] 

In the next Section we examine the vortex/anti- vortex background and argue that the 
two-state, two-dimensional, fermion parity preserving realization can be established. The 
physical picture that emerges is that there are two towers of states, one built on an "empty" 
zero energy state \b), the other on the "filled" zero-energy state [/), and the A operator, 
which connects the two "vacua," fills or empties the zero energy state. 




(44) 
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Vort ex/ant i- Vortex Order Parameter 



Insight on physical states in the presence of a vortex in a superconductor adjoined to a 
topological insulator can be gotten by considering a vortex/anti- vortex background. The 
zero energy mode for an anti-vortex at the origin, A(r) = v(r)e~ ld , is given by 



= N J (/ir)exp{ivr/4- F(r)} 
^ = NJi (//r)exp{-i(# + 7r/4) -V(r)} 



(45) 



To simplify the discussion, we omit the chemical potential and evaluate V(r) = f r dr' v(r') 
with the asymptotic form of v(r) m. Thus the zero-energy mode for the vortex 

becomes, approximately 

% ^Ne-™' 4 e- mr ( M (46a) 



(46b) 



V J 

while the anti-vortex at r = R, in the same approximation leads to 



-m|r-R| 



The corresponding 4-spinors that solve Eq. (11) at zero energy are 

\ 



( Ne-^/^e 



V —N e* 71 "/ 4 e~ mr j 
( \ 

e -i7r/4 g -m|r-R| 




(47a) 



(47b) 



Consider now a configuration with a vortex at the origin and an anti- vortex at R. No 
zero mode is present in the spectrum of h; rather there are two bound states, one with 
positive, exponentially small energy e e -mR anc j trig other with equal magnitude, but 
opposite sign. 

The former, called $" v , consists of portions localized at the origin (vortex) and at 
r = R (anti- vortex) . The latter is given by <3?™ = C$™*, and has similar structure. Both 
contribute unambiguously to the expansion of the quantum field operator Vl/, the former 
with an annihilation operator, the latter with a creation operator. 



* = ^cont + ae e~ i£t $f + a\ e iet C<S>\ 



(48) 
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The first term on the right is the continuum contribution, as in (40). The Fock space 
spectrum is clear. There is a vacuum state |0) annihilated by a e 



a e \n} = 0. (49a) 
A low- lying state is gotten by operating on |f2) with a e . 

4 |n) = |/) (49b) 

a e |/) = |fi) (49c) 
a\\f) = (49d) 

The remaining states, created by a E can be built either on the vacuum |n) : a E a E/ a E// ... |n), 
or on the low lying state |/) : a\ |n): a E a E/ a El/ . . . \f). 

Now let us remove the anti- vortex by passing R to infinity. Both \P!j? e collapse to their 
zero-mode limit, ty± e Vl/p, and the expansion (42) becomes 

* = *cont + (^f) V2^ 

V V2 ' (50) 

= *cont+^V^*g 

Moreover the action of A = -k= (a e + a e ) may be read off (49). Renaming |n) as \b), 
we find 

{A, A} = 1 

A\b) = ± |/) (51) 
A |/) = ^ |6) 

and two towers of states are built upon | b) and | /) . 

In this way we justify the two-dimensional, fermion parity preserving realization of the 
zero mode algebra in a superconducting/topological insulator system. 

[Note the occurrence of the factor y/2 modifying ^q- This explains its first appearance 
in eq. (40). This factor compensates in the completeness sum for the loss of the anti- vortex 
wave function.] 

Because no explicit solutions in a vortex/anti- vortex background are available, the 
argument in this Section is qualitative, without explicit formulas. However, one may 
consider a one-dimensional example with Major ana fermions in the presence of a kink 
and/or a kink anti-kink pair. [21] In that model one can solve equations explicitly and 
verify the behavior described here for the two-dimensional vortex case. In this way one 
also establishes that even in one spatial dimension (in the absence of rotation and spin 
to enforce fermion parity) the two-dimensional realization of the zero mode algebra is 
appropriate. 
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In Appendix B we present an approximate determination of the low-energy eigenvalues 
in the presence of a vortex/anti- vortex pair. The result supports the above qualitative 
argument: an exponentially small splitting of the zero-energy mode is established. Also in 
the Appendix, we study the two vortex background, and find, within the same approxima- 
tion that no energy splitting occurs; rather two zero modes persist as anticipated by index 
theorems. 
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Quantizing Majoranas Fermions in Generic Superconductors 



In the preceding sections we showed that Majorana's quantization prescription of the Dirac 
equation directly applies to the the full quantum field describing the proximity effect of an 
s-wave superconductor to surface states of a topological insulator. Below we shall show that 
Majorana's quantization prescription of real neutral fermions is rather generic in super- 
conductors, with or without Dirac-type dispersions. The construction below is possible for 
any half- integer spin (fermionic) particle. The reality conditions on the fermionic fields fol- 
low from symmetries of the Bogoliubov-de Gennes (BdG) Hamiltonian for superconductors 
constructed in the Nambu basis. 

Let us consider a system with fermionic degrees of freedom if) r ns an d V4,n,s> where r 
labels position, n the possible flavors (bands, for instance), and s the spin (half-integer) 
along a chosen quantization axis. For simplicity, we shall define an index a = (r, n, s) 
that encodes all these degrees of freedom. The Hamiltonian describing superconductivity 
in such a system can be written as 

a,/3 

= E \^ H *P ^ + {~H T ) ap i>l + -iPl A aP i>l + ^Va V/3 • (52) 



a,/3 



Defining 



we can write 



n = \ ¥ ( ft _% * ) * = & h * ■ ( 54 ) 

That H T = H* follows from H = W . Notice that A = — A T is enforced because of 
fermionic statistics, and consequently we can also write 

H-\(_l. _%)■ (55) 

Let us define 

C=(° T i) , (56) 



x 1 , 

so that C = C T = C* = = C~ l . The operators \l/ must satisfy the constraint 

Cab *J = *a , (57) 
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where the index a = (a,p), with p = ± the Nambu grading (^ a ,+ = ipa and ^ a ,- = Vi)- 
The fermionic commutation relations of the fields ip, ip' translate into 



{# a ,* 6 } = C a6 and {^ a M} = S ab . (58) 



Conjugation symmetry 

One can easily check that any BdG-type h as in Eq. (55) possesses the following conjugation 
symmetry, 

-h*=C*hC. (59) 

We stress that fermionic statistics underlies this result, as it is the reason for the minus 
signs and the complex conjugation in both terms in the second row of Eq. (55). 

It follows from this symmetry that positive and negative eigen modes of h are paired: 

h<5> E = E$ E h(C$* E ) = -E(C<f>* E ) , (60) 

or equivalently 

c$; E = $_ B . (6i) 



Generic Majorana basis and its real equation of motion 

Consider a unitary transformation V, under which 

h^h = VhV^. (62) 

It follows that 

-h* = V* {-h*)V* ] = V* {C*hC)V* ] 

= V* C*V ] VhV ] VCV* ] 

= (VCV T )* h (VCV T ) 

= C*hC, (63) 

so the transformation law of C is 

C ^ C = VCV T . (64) 

(Notice that CC* = VCV T V* C* W = /, so C' 1 = C* still.) 

We will construct below a unitary matrix V such that C = I. This basis is the 
appropriate Majorana representation for the generic superconducting system of half-integer 
spin particles (for any number of flavors). In this basis, one has h = —h*, so that h is 
imaginary, or equivalently ih is real. It follows from Schrodinger's equation that 



dt + ih)V = 0, (65) 
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so the equation of motion for the field is purely real and thus admits purely real solutions. 
Notice that this path mirrors Majorana's formulation of the Dirac equation for spin 1/2 
particles (he constructed a purely imaginary representation of the Dirac matrices, obtaining 
an equation of motion that was real). 

Notice that in this basis the commutation relations become 

{^ a , * h } = C ab = 5 ab and {§ a , *£} = 5 ab , (66) 

corresponding to real fermions 

*a = • (67) 
Construction of V 

The appropriate unitary transformation V which makes C = I is constructed as follows. 
Because of fermionic statistics, the time-reversal operator O squares to —1. One can write 
= TIC, where K, is complex conjugation and T = e mSV , with S y the y-component of 
the angular momentum operator (in a representation such that S y is a purely imaginary 
matrix). T is a real anti-symmetric matrix (7~ = T* and = —7"), with T 2 = e l27TSV = 
— -^spinxflavor when spin is half-integer. For instance, for spin 1/2 particles T = ia 2 . 
Consider the following transformation. 

v= U; "q- + ) £ " /4 ' <*±A<- ijftr > (68) 

[compare with eq. (26).] Notice that Q± are projectors {Q"± = Q±), that Qj\_ = Q±, and 
that Q\ + QPl = Ispinxflavor and Q + Q- = Q~Q+ = 0. Also notice that Q*± = Q± = Q T . 
One can then easily check that the above defined V is such that 

C = VCV T = 1 . (69) 
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Summary 



In this paper we studied mainly three issues regarding the quantization of Majorana 
fermions in superconductors, following closely Majorana's original definitions, and looked 
beyond just the Majorana zero energy modes that are bound to topological defects such 
as vortices. 

We started by analyzing the specific case of Dirac-type systems describing s-wave in- 
duced superconductivity on the surface of topological insulators. We showed that the entire 
ijj field of the superconductor model (and not merely particular modes) obeys equations 
that are analogous to the Majorana equations of particle physics. 

We then analyzed the quantization of the theory in the presence of vortices. We showed 
that fermion parity can be preserved, even with a single zero energy state. This quanti- 
zation scheme shows that one can obtain a two-dimensional Hilbert in the presence of a 
single vortex in an infinite plane, presenting a case where each Majorana fermions can be, 
when present in odd numbers, more than "half" a qubit. 

Finally, we showed that the Majorana quantization procedure that we discussed for 
the Dirac-type equations describing s-wave induced superconductivity on the surface of 
topological insulators does extend, more broadly, to any superconductor. The constraints 
imposed by fermionic statistics on the symmetries of Bogoliubov-de Gennes Hamiltonians 
are sufficient to allow real field solutions in the constrained doubled Nambu space that can 
then be quantized as Majorana fields. This results follows simply from fermionic statistics 
plus superconductivity, irrespectively of the presence or absence of any other symmetries 
in the problem, such as spin rotation invariance or time-reversal symmetry. 
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Appendix A 

We present the 4-component, positive energy solutions to (11). The eigenvalues 



E± = vWmF + ^F 

are associated with the eigenvectors 



$ + (k)- ' 



^E + -k+^ |+) 
A i \ 



$ (k) _ 1 I 



The negative energy spinors are given by C$^_ (— k). |±) are defined in the text. 

Appendix B 

We study the low-lying energy levels of the Dirac-type Hamiltonian h in (9) with /u set to 
zero and A chosen first in an approximate vortex/anti- vortex profile, 

A„ c = m e^ r ~ R / 2 ) e -^(r+R/2) (R1) 
and then similarly with two vortices. 

A OT = me ffl ( r - R «e ffl ( r+a / 2 ' (B.2) 
Here O is the argument of the appropriate vector, 

e iQ(r) = x+iy = e i9 

e i«(±R) = ± (wy) _ _j_ e ie 

with x = rcos6,y = rs'm9,X = RcosO,Y = i?sin0. One vertex is located at r ~ R/2, 
the anti- vortex or the second vortex at r ~ — R/2. [23] 

B.l Vortex/anti- Vortex 

Near the vortex at r ~ R/2 the order parameter A vi) is approximated by 

A„ e ->A„ = me in(r-R/2) e -iQ(R) 
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The zero mode in the presence of A„ differs from (47a) by a phase, due to the additional 
phase e t& in A v . Also the location is shifted by R/2. 



to 



( v \ 






m 

/7 



, -i(vr/4+e/2) e -m|r-R/2| 



(B.4) 



Similarly, with the order parameter near the anti- vortex at r ~ —R/2 taken as 

A e „->A e = me in(-R) e -in(r+R/ 2 ) 

= _ m e iB e -iU(r+-R/2) 
the zero mode solution replacing (47b) reads 

/ \ 



(B.5) 



r 



V 

v* 

V o / 



V = — e^T+Q/ 2 ) e -m|r+R/2| 



7T 



(B.6) 



Next we evaluate the matrix element between ipQ and Vo of h in (9), with fx = 
and order parameters as in (B.l). We find that the diagonal matrix elements vanish 
{v\ h \v) = (v\ h \v) = 0. The energy shift A E, determined by the off-diagonal elements of 
h, is 



AE = ±| (v\h\v) 
{v\ h \v) = (v\ h \v)* = 



d 2 r (v*p + v + v A* s v) — h.c. 



(B.7) 
(B.8) 



The evaluation of the first integrand proceeds by recalling that p+v = A v (r — R/2) v* 
and yields, after a shift of integration variable by R/2, 



d 2 rv*p + v = -iE*r / rdre~ mr \ d6 e i( - e ~ e ^ 
n Jo J-ir 
exp - m[r 2 + R 2 + 2rRcos(6 - 9)] 1/2 . 

A further shift of 9 by leaves 



(B.9a) 



d rv* p + v 



■2m 3 



r dr e 



dO cos 9 e- mD 



(B.9b) 



D = Vr 2 + R 2 + 2Rr cosO 
For the second integrand, a similar shift, first by R/2 and then by gives 



r dr e 



r + RcosO D 
d9 e 



(B.lOa) 
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We notice that the 9 integrand may also be presented as — ^ e mD , thereby transforming 
(B.lOa) after an integration by parts into 

r pOO p7T 

d 2 rvA* V vV = i^ / dr{\ - rm)e~ mr d6 e~ mD . (B.lOb) 
J Jo Jo 



Thus we find that 

(v\ h \v) = ie 



e = Am 2 i dre- mr \\ dO[l - rm{\ + cos 9)} e~ mD 
Jo 71 Jo 

and the energy is shifted from zero by ± e. 

Numerical integration of (B.ll) at large R yields a result consistent with 



(B.ll) 



H - m (mR) 1/2 e~ mR . (B.12) 



R-^OO y ^ 

This may be derived analytically with the following argument. We replace the upper limit 
(oo) of the r integral by R and approximate D by R + r cos 6. The 9 integral now leads to 
modified Bessel functions Jo and I±, and we keep only their large argument, exponential 
asymptote. The remaining r integral yields (B.12). 

B.2 Two Vortices 

The order parameter (B.2) describing two vortices located at r ~ ± R/2 reduces at r ~ R/2 
to 

A VV ^A V+ = me ai(r-R/2) e in(R) 

(B.13a) 

while the one at r ~ —R/2 becomes 

A ra ^A,_ = me in(-R) e in(r+R/2) 

= _ m e ie e if7(r+R/2) _ (B.13b) 

The corresponding zero modes differ by phases from the vortex solution (47a) or (B.4) but 
they retain their spinor structure. 



r 



( v + \ 




V -< / 
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r 



(B.14) 



V 



The explicit expressions for u + and V- are not needed, because the above form of the 
spinors guarantees that all matrix elements of h vanish. Thus, within our approxima- 
tion, the 2-vortex background retains its two zero modes. This is to be expected because 
asymptotically such a configuration is indistinguishable from a double vortex, 



A 



vv r+oa 



me in(r) e iO(r) 



m e 



2 it 



(B.15) 



and a double vortex possesses two zero modes. [1] 
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